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Abstract 

Chiral rotation observed in 128 Cs is studied using the newly developed 
microscopic triaxial projected shell model (TPSM) approach. The observed 
energy levels and the electromagnetic transition probabilities of the nearly 
degenerate chiral dipole bands in this isotope are well reproduced by the 
present model. This demonstrates the broad applicability of the TPSM ap- 
proach, based on a schematic interaction and angular-momentum projection 
technique, to explain a variety of low- and high-spin phenomena in triaxial 
rotating nuclei. 

The classification of band structures from symmetry considerations has 
played a central role in our understanding of nuclear structure physics. Most 
of the rotational nuclei are axially symmetric with conserved angular-momentum 
projection along the symmetry axis. This symmetry has allowed to classify a 
multitude of rotational bands using Nilsson scheme and has been instrumen- 
tal to unravel the intrinsic structures of deformed nuclei [H, 0]. Although, 
most of the deformed nuclei obey axial symmetry at low-excitation ener- 
gies and spin, there are also known regions of the periodic table, referred 
to as transitional nuclei, that violate axial symmetry and are described us- 
ing the triaxial deformed mean-field. Further, some nuclei that are axial 
in the ground-state become triaxial at higher excitation energies and angu- 
lar momenta. There are several empirical observations indicating that axial 
symmetry is broken in transitional regions. For instance, the moment of 
inertia of the transitional nuclei rises sharply with angular-momentum for 
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the ground-state band that cannot be explained using axial symmetry [3j. 
Further, the transition quadrupole moment of the ground-state band of tran- 
sitional nuclei increases with spin and has been explained using the triaxial 
mean- field model Q. We would like to mention that this is still a debateable 
issue, these features could also be explained using the 7-vibrational degree of 
freedom rather than the stable 7-deformation used in the above cited work. 

Recently, it has been shown that rotational motion of a triaxial deformed 
mean-field can also exhibit chiral symmetry with angular-momenta of va- 
lence protons, valence neutrons and the triaxial core directed along the three 
principle axes of the triaxial field. Chiral symmetry emerges from the combi- 
nation of triaxial geometry with an axis of rotation that lies out of the three 
principle planes of the ellipsoid [5] . It has been shown that nuclei in the mass 
~ 130 region exhibit chiral geometry with proton in the lower-half of the /in/2 
subshell with its angular-momentum aligned along the short-axis, neutron in 
the upper-half of the subshell with angular- momentum along long-axis and 
the core angular-momentum directed along the medium axis as the energy 
is minimum for the triaxial shape along this axis 0, 0, 0, 0, lioj - Till date, 
the best candidates for chiral doublet band structures are observed in 128 Cs 
and 135 Nd isotopes ll|, ll3j. The measurement of the lifetime of the states 



in the degenerate dipole bands plays an extremely crucial role in identifying 



the chiral rotations in atomic nuclei [14| and these measurements have been 
performed for 128 Cs. 

Several theoretical approaches have been used toprobe the chiral doublet 
band structures observed in odd-odd nuclei 15, 16, 17, 18, 19, 2(|. In most of 



the analysis, particle-rotor approach with valence proton and neutron coupled 
to the triaxial rotor has been employed. Although, this model has been shown 
to provide a reasonable description of the most of the characteristics of the 



chiral bands, it is known to depict discrepancies with the observed data [11 
For instance, the measured B(E2) values for the yrast band in 128 Cs drop with 
spin and the two-quasiparticle triaxial rotor (2QPTR) calculations display 
opposite trend of increasing B(E2) with spin [ll]. The discrepancy is possibly 
related to the phenomenological nature of the particle-rotor model approach 
with the triaxial core held fixed for all spin values. In the present TPSM 
approach, the core or vacuum state has a microscopic structure as compared 
to 2QPTR where it is simply a triaxial rotor. 

The purpose of the present work is to provide a detailed investigation 
of the chiral rotation using the microscopic triaxial projected shell model 
(TPSM) approach. This model has been recently used to investigate 7- 
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vibrational and quasiparticle bands in even-even and odd-mass triaxial nu- 
clei and it has been demonstrated to provide an accurate description of the 



observed properties |2l|, |22|, |23|, |25|, |26|. In the present study, we have 



generalized the TPSM approach for odd-odd mass and as a first application 
of this development, the chiral rotation observed in 128 Cs is investigated. 
In particular, the transition probabilities, which contain the important in- 
formation on the chiral geometry, are evaluated in the present work. The 
angular momentum projection technique used in the present model provides 
the B(E2) and B(M1) transition probabilities between states with well de- 
fined angular momenta and a direct comparison with the measured values is 
feasible. 

The basis space of the TPSM approach for odd-odd nuclei, developed in 
the present work, is composed of one-neutron and one-proton quasiparticle 
configurations: 

{\<f> K >=a v W\0>}. (1) 

The above basis space is adequate to describe the ground-state configuration 
of odd-odd nuclei. For the investigation of the excited configurations, proton- 
and neutron- pairs need to be broken, which shall form focus of our inves- 
tigation in future. The triaxial quasi-particle (qp)-vacuum |0 > in Eq. §T§ 
is determined by diagonalization of the deformed Nilsson Hamiltonian and a 
subsequent BCS calculations. This defines the Nilsson+BCS triaxial qp-basis 
in the present work. The number of basis configurations depend on the num- 
ber of levels near the respective Fermi levels of protons and neutrons. The 
configuration space is obviously large in this compared to the nearby 

odd-mass nuclei, and usually several configurations contribute to the shell 
model wave function of a state with nearly equal weightage. This makes 
the numerical results very sensitive to the shell filling and the theoretical 
predictions for doubly-odd nuclei become more challenging. 

The states \<p K > obtained from the deformed Nilsson calculations do not 
conserve rotational symmetry. To restore this symmetry, angular-momentum 
projection technique is applied. The effect of rotation is totally described 
by the angular-momentum projection operator and the whole dependence 
of wave functions on spin is contained in the eigenvectors, since the Nilsson 
quasi-particle basis is not spin-dependent. From each intrinsic state, k, in (pQ) 
a band is generated through projection technique. The interaction between 
different bands with a given spin is taken into account by diagonalising the 
shell model Hamiltonian in the projected basis. 
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The Hamiltonian used in the present work is 

H = H - l - X QlQv - G mP ] P ~GqJ2 P IK ( 2 ) 
and the corresponding triaxial Nilsson Hamiltonian is given by 

Hn = h - | e g + A 2 ^- 2 j , (3) 

where Ho is the spherical sing le-particle shell model Hamiltonian, which con- 



tains the spin-orbit force [27| . The second, third and fourth terms in Eq. ([2]) 
represent quadrupole-quadrupole, monopole-pairing, and quadrupole-pairing 
interactions, respectively. The axial and triaxial terms of the Nilsson poten- 
tial in Eq. [3] contain the parameters e and e', respectively, which are re- 
lated to the 7-deformation parameter by 7 = tan" 1 — . The strength of the 
quadrupole-quadrupole force x is determined in such a way that the em- 
ployed quadrupole deformation e is same as obtained by the HFB procedure. 
The monopole-pairing force constants Gm used in the calculations are 

N — Z 

G V M = [20.12 - 13.13 - ]A-\ G n M = 20.12A- 1 . (4) 

Finally, the quadrupole pairing strength Gq is assumed to be proportional 
to the monopole strength, Gq = 0A6Gm- All these interaction strengths are 
the same as those used in our earlier studies for the even-even nuclei in the 



same mass region |3|, |22|, |23j, |24j. Thus, we have a consistent description for 
doubly-even and doubly-odd nuclear systems. 

Once the projected basis is prepared, we diagonalize the Hamiltonian in 
the shell model space spanned by Pm K \4>k > ■ The projected TPSM wave 
function is given by 

\a,IM> = Y<K pI MK\<i>K> ■ (5) 

K,K 

Here, the index a labels the states with same angular momentum and k 
the basis states. In Eq. f° are the weights of the basis state k, P[ ik is 



three-dimensional angular momentum projection operator [28 

Pmk = J dn D J MK (n) R(Q). (6) 
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Finally, the minimization of the projected energy with respect to the expan- 
sion coefficient, f%, leads to the Hill- Wheeler type equation 

J2{H KK >-E a N KK ,)f° = 0, (7) 

K> 

where the normalization is chosen such that 

Y^r.N^rj = 5 aa , (8) 

KK-' 

The angular-momentum-projected wave functions are in the laboratory 
frame of reference and can thus be directly used to compute the observables. 
The reduced electric quadrupole transition probability B(E2) from an initial 



state (o"j, Ij) to a final state (cr/, If) is given by [29 



e 2 



B{E2J l ^I f ) = w ^- l \ <t7 / ,/ / ||Q a ||t7 i ,/ i > | 2 . (9) 

In the calculations, we have used the effective charges of 1.6e for protons 
and 0.6e for neutrons. The reduced magnetic dipole transition probability 
B(M1) is computed by 

B(Ml,I r ^I f ) = ^ 1 \ «T f ,I f \\Mi\\<Ti,Ii> | 2 , (10) 
where the magnetic dipole operator is defined as 

Ml = gff + (91 - gJW- (ii) 



Here, r is either v or ir, and g\ and g s are the orbital and the spin gy- 
romagnetic factors, respectively. In the calculations we use for gi the free 
values and for g s the free values damped by a 0.85 factor, i.e., 

g? = l, g? = 0, ^ = 5.586 x 0.85, ^ = -3.826 x 0.85. (12) 

Since the configuration space is large enough to develop collectivity, we do 
not use core contribution. The reduced matrix element of an operator O (O 
is either Q or M) is expressed by 
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Figure 1: (Color online) The angular- momentum projected bands obtained for different 
intrinsic K-configuration, given in legend box, for 128 Cs nucleus. 



In the present study, a detailed investigation of the chiral band structures 
observed in 128 Cs has been performed using the TPSM approach. We have 
chosen this nucleus since it is the only system for which exhaustive lifetime 
measurements have been performed. As already discussed, the electromag- 
netic transition probabilities are important to probe the chiral symmetry 
breaking of the observed dipole bands. In Fig. 1, the projected bands, ob- 
tained from the triaxially deformed intrinsic Nilsson state by performing the 



<(Tf,If\\0 L \\(Ti,Ii> 



J IiKiJ 1 



.'°7 
I f K f 




(13) 



6 




Figure 2: (Color online) Comparison of the TPSM energies after configuration mixing with 
the available experimental data for the yrast and side bands of the studied 128 Cs nucleus 
is shown in panel (a). Calculated staggering parameter S(I) = [E(I) — E(I — 1)] /2J is 
plotted along with the measured values in panel (b). 



three-dimensional angular-momentum projection, are displayed. Nilsson in- 
trinsic states have been obtained with deformation parameters, e = 0.220 and 
e' = 0.14. The axial deformation value has been adopted from particle- rotor 
model and other studies and the chosen value of non-axial deformation corre- 
sponds to 7 ~ 30°. The lowest bands correspond to the projected states from 
the intrinsic configuration with quasiparticle energy of 1.34 MeV and are ob- 
tained by specifying K-values in the three-dimensional angular-momentum 
projection operator, given in Eq. (EJ). The projected bands obtained with 
K=4, 5 and 7 are degenerate for 1=9 and for higher spins up to 1=13, K=5 
and 7 are degenerate with K=4 band becoming less favoured with increasing 
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•— « e' = 0.0250 
-< — < e = 0.0750 




Spin (II) 



Figure 3: (Color online) Calculated energy difference (AE) between the yrast and side 
bands as a function of angular-momentum for a range of the triaxial deformation e' values. 



spin. For 1=14, it is noted that bands with K=l and 2 belonging to the 
quasiparticle configuration of energy 1.77 MeV cross the bands belonging to 
the quasiparticle energy of 1.34 MeV. It is further observed that more bands 
belonging to the higher quasiparticle energy cross at higher spins. This cross- 
ing of the quasiparticle states has not been explored in the earlier studies of 
the chiral bands. 

In the second stage of the calculations, projected bands obtained above 
are used to diagonalize the shell model Hamiltonian, Eq. (j2]). It is important 
to mention that the projected bands are not orthogonal, in general, and 
need to be orthogonalized before the diagonalizatoion of the Hamiltonian. 
The method that is adopted in the present approach, and in most of the 
generator coordinate methods, is to calculate the eigen-values of the norm 
matrix. The norm eigenvalues are non-negative quantities since the norm is 
a positive semi-definite matrix. However, it is quite possible that some of 
them may vanish. This happens under the circumstances that the multi-qp 
states become linearly dependent (non-orthogonal) when they are projected. 
Such redundant states are removed simply by discarding the zero eigenvalue 
solutions of the norm eigen-value equation. The details are given in the 
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appendix of ref. [30 



Bands after diagonalization have dominant compenents from several angular- 
momentum projected basis configurations and the lowest two bands obtained 
are displayed in the top panel of Fig. 2 along with the corresponding experi- 
mental chiral bands known for 128 Cs. It is evident from Fig. 2(a) that TPSM 
calculations reproduce the known experimental energies of the chiral bands 
quite well. It is noted from the figure that the two bands are very close to each 
other up to 1=14, and after this spin the two bands deviate with increasing 
spin. In Fig. 2(b), the staggering parameter, S(I) = [E(I) — E(I — 1)] /2J 
for the two calculated chiral bands are plotted and compared with the corre- 
sponding experimental numbers. It is noted that for low angular-momenta, 
the calculated quantity deviates considerably from the experimental val- 
ues, but for higher angular-momenta a reasonable agreement is obtained. 
Particle-rotor model results depict a similar discrepancy at lower angular- 
momenta 0, [l5| . 

To elucidate the importance of the triaxial deformation in the occurence 
of degenerate dipole bands, in Fig. 3 the energy difference (AE) between the 
two dipole bands is drawn as a function of the the angular-momentum for a 
range of the triaxial deformation values. It is evident from the figure that 
for lower values of e', close to the axial limit, AE is quite large, especially 
for higher angular-momentum states. However, for e' larger than 0.125, it is 
noted that AE becomes quite small for all angular-momentum states. 

We now discuss the electromagnetic transition probabilities of the ob- 
served doublet bands in 128 Cs. B(E2) (from I to 1-2) and B(M1) (from I to 
1-1) transition probabilities have been calculated using the expressions, Eqs. 
(Q and fflOl) and are shown in Fig. 4 for both the yrast and the side bands. 
B(E2) transition probabilities for the yrast band, displayed in Fig. 4(a), de- 
pict a slight drop at around 1=13 and is due to the crossing of the quasiparti- 
cle state with energy equal to 1.77 MeV with the ground-state quasiparticle 
configuration. The measured B(E2) for the yrast band, displayed in Fig. 4(a), 
shows a more pronounced drop at a slightly higher spin value of 1=15. Fur- 
ther, the measured B(E2) after 1=16 show a smooth behaviour and is well 
reproduced by the TPSM calculations (see Table I). The calculated B(E2) 
for the side band displays a little odd-even staggering, in particular, for in- 
termediate spin values. Except for 1=13, the measured B(E2) for the few 
transitions are in good agreement with the calculated values and is quite 
evident from Table I. 

The present TPSM approach provides a little improved description of the 
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Table 1: Experimental and calculated intra band B(E2) transition strengths (from 
I I — 2) for yrast and side bands of 128 Cs in (W.u.) 



1(h) 


B(E2) 
Yrast Band 




B(E2) 
Side Band 






Experiment 


Theory 


Experiment 


Theory 


11 




25.19 




25.19 


12 




26.10 




22.37 


13 


48+21 


22.94 


qo+63 


32.46 


14 


73+}? 


26.29 




19.84 


15 


35±S° 


25.90 


42 + g 4 


39.49 


16 


42+* 2 


26.90 


13 5 +6 -° 


14.18 


17 


36+* 6 


34.47 


34.58 


18 


35.0^° 


32.14 




28.08 


19 


30 + s 


35.10 




22.57 


20 


2011 


22.40 




33.02 



measured B(E2) values due to its microscopic treatment of the core (or vac- 
uum) as compared to the phenomenological rotor core in 2QPTR. In TPSM, 
vacuum configuration (quasiparticle states) are generated by first solving tri- 
axial Nilsson potential with the parameters e and e'. In the present study, 
we have adopted e = 0.22, which is, as a matter of fact, same as that em- 
ployed in the particle-rotor model analysis, e' = 0.14 has been used as this 
corresponds to 7 ~ 30° and is important to lead to chiral configuration. In 
the second phase, standard BCS calculations are performed to obtain the 
quasiparticle states with the Nilsson basis. Although, in TPSM angular- 
momentum projection is performed with pairing (monopole + quadrupole) 
plus quadrupole-quadrupole interaction for each angular-momentum, but due 
to limited quasiparticle space with only one-proton and one-neutron excita- 
tions considered, it is expected that deformation of the projected states will 
not be very different from the vacuum or mean-field configuration. In TPSM, 
the vacuum configuration or mean- field is held fixed for all the projected 
states and this is certainly a drawback of the present analysis. In a more 
accurate treatment, projection after variation, mean-field is re-calculated for 
each angular-momentum state. We consider that due to this problem in the 
present approach, we observe only a marginal drop in the B(E2) in compar- 
ison to the experimental data. 

B(M1), shown in Fig. 4(c), for the yrast band are initially smooth, but 
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after spin, 1=13 show large odd-even staggering. For the side band, B(M1) 
has slightly a different behavior with B(M1) transitions showing a drop for 
1=13 to 16. The corresponding measured transitions are almost constant for 
1=13,14 and 15, and depict a drop for 1=16. The phase of the odd-even 
staggering and the magnitudes of the B(M1) transitions are quite similar in 
the two dipole bands and is one of the important criterion for the chiralilty 
of the two bands. The calculated B(M1)/B(E2) ratios, shown in Figs. (4e) 
and (4f), are in good agreement with the experimental values obtained from 
the intensities [7J. The inter-band B(M1) transitions between yrast to side 
and side to yrast, shown in Figs. 4(g) and 4(h), have same phase of odd- 
even staggering and further are similar in magnitude. What is important is 
that this phase is opposite to that of the intra-band B(M1) transitions of 
Fig. 4(c) and 4(d) and forms another essential benchmark for the existence 
of the chiral geometry j^j. In Fig. 4(h), the measured B(M1) values for a few 
available transitions are noted to be in good agreement with the calculated 
values. 

In conclusion, the chiral dipole band structures observed in 128 Cs have 
been investigated using the triaxial projected shell model for odd-odd nuclei. 
It has been demonstrated that observed energy levels and electromagnetic 
transition probabilities of the dipole bands are quite well reproduced. Fur- 
thermore, it has been demonstrated that the calculated inter band B(M1) 
transitions from yrast to side and vice-versa, depict opposite staggering be- 
haviour as compared to the the in-band transitions and are in conformity 
with the predicted properties of the chiral bands. 

The authors would like to acknowledge Dr. Ernest Grodner for providing 
the measured data on the electromagnetic transitions. 
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Figure 4: (Color online) Comparison of the experimental B(E2), B(M1) transition 
strengths and their ratio for 128 Cs [Txj, [l2| with the calculated values for yrast and side 
bands are depicted in panels (a) to (f). The calculated interband transitions are shown in 
panels (g) and (h), and panal (h) also depicts the measured values for a few transitions. 
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